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Abstract 

The influence of the dimerization process on the nematic ordering is 
investigated by using a nonextensive thermostatistics, namely Tsallis thermo- 
statistics(TT). A theoretical model taking into account the dimerization in- 
fluence on the nematic scalar order parameter has been recently presented[l]. 
Our study considers this model and the theoretical predictions for the nematic 
order parameter are improved by using TT. 
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I. INTRODUCTION 



Since the papers by Tsallis [2,3], Tsallis thermostatistics (TT) has been commonly 
used for the investigation of the physical systems. Boltzmann-Gibbs statistics is a powerfuU 
tool to study a variety of the physical systems. However it fails for the systems which i) have 
long-range interactions, ii) have long-range memory effects, and in) evolve in multi-fractal 
space-time. The system which has the above properties is called " nonextensive one". If 
investigating a nonextensive system, one then needs an appropriate statistics, and TT is one 
of these statistics. It is the aim of this paper to enlarge the apphcations of TT to the hquid 
crystal systems. In doing so, we first summarize the axioms of TT, and then investigate the 
influence of the dimerization process on the nematic ordering by using TT. 

TT has been apphed to various concepts of thermostatistics and achieved in solving 
some physical systems, where Boltzmann-Gibbs statistics is known to fail [4]. Recently, 
MST has been generalized within TT and generalized MST has been applied to a nematic 
liquid crystal, p-azoxyanisole (PAA), in [5] in which Kayacan et al. used the second choice 
for the internal energy constraint which will be given below. 

TT considers three possible choices for the form of a nonextensive expectation value. 
These choices have been studied in [6] and apphed to two systems; the classical harmonic 
oscillator and the quantum harmonic oscillator. In that study, TsaUis et al. studied three 
different alternatives for the internal energy constraint. The first choice is the conventional 
one and used in [2] by TsaUis 
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l:P^ei^U('\ (1) 

i=l 

The second choice is given by 

w 

T.i^e. = UlP (2) 

i=l 

and regarded as the canonical one. Both of these choices have been applied to many different 
systems in the last years [7]. However both of them have undesirable difficulties. The third 
choice for the internal energy constraint is 
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This choice is commonly considered to study physical systems because it is the most appro- 
priate one, and is denoted as the Tsallis-Mendes-Plastino (TMP) choice, q index is called 
the entropic index and comes from the entropy definition [2] , 

Sq = k ' (4) 

where A; is a constant, X^iPi = 1 is the probability of the system in the i microstate, W is the 
total number of configurations. In the limit g ^ 1, the entropy reduces to the well-known 
Boltzmann-Gibbs (Shannon) entropy. 

The optimization of Sq leads to 
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This equation is an implicit one for the probabilities pj. Therefore the normalized q — 
expectation value of an observable is defined as 

A, = = (A.), (7) 

Z^j=l Pi 

where A denotes any observable quantity which commutes with the Hamiltonian. This 
expectation value recovers the conventional expectation one when g = 1. As mentioned 
above, Eq.(7) is an implicit one and in order to solve this equation, Tsallis et al. suggested 
two different approaches; "iterative procedure" and "/9 — > p'" transformation. 

Since the present study is based on a theoretical model [1], which has been recently 
presented, at this stage we give a summary of this model. 

The contituent molecule in the liquid crystals is assumed as a rigid rod in the most 
commonly used theoretical models. Maicr-Saupe theory [8,9] and this approximation ex- 
plain the most of the macroscopic observable in nematics. However, this approximation not 



appropriate to explain the observed anomalies in 4,n-alkyloxybenzoic acids [10-12]. They 
have nematic and smectic phases. The anomalies are in the temperature dependence of the 
elastic constants[13], the electroconductivity [14], and the dielectric permitivity [15]. The 
variation of these constants strongly deviate from those of the classical nematics. The un- 
usual thermal behaviour of the nematic phase of 4,n-alkyloxybenzoic acids was found and 
investigated by some optical methods. The experimental results show a strong textural tran- 
sition and an increasing of the depolarized Rayleigh scattering at a given temperature T* 
in the nematic range [11]. It is reasonable to assume that this phenomenon is connected to 
the temperature variable form of the constituent molecules of the 4,n-alkyloxybenzoic acids. 
These substances change the molecular structure with the temperature between three states, 
which are closed dimers, open dimers and monomers [16,17]. The variations of the closed 
dimer, open dimer and monomer equilibrium concentrations vs temperature were detected 
by FT-IR spectroscopy [17]. On the dimer ring, realized by the pairing of two monomers 
in a closed dimer by hydrogen bonding, are imposed different conformational states, which 
were detected by far FT-IR spectroscopy [18-19]. A deviation of the temperature trend 
of the order parameter in the nematic phase of the NOBA (nonyloxybenzoic acids) was 
demonstrated [20]. A similar temperature trend for NOBA was found in [21]. All these 
experimental results indicate that the dimarization process strongly influences the nematic 
ordering in the 4,n-alkyloxybenzoic acids. 

Now, we discuss the classical theories. A model keeping the mean field approximation, 
but considering the liquid crystal molecules as pairs of monomers [22] or as association of 
structural elements in flexible molecules with internal degrees of freedom [23] were presented. 
These theories are motivated by the insufficiency of the classical models, which approximate 
the liquid crystal molecules as rigid rods, to explain the phenomena like reentrant liquid 
crystal state modeled [22] and the isostructural phase transitions [23]. Barbero et al. found 
a similarity in the liquid crystal state modeled [22] as a mixture of monomers and dimers, 
used to explain the provocation of the reentrant behaviour and the nematic state of the 
4,n-alkyloxybenzoic acids detected by FT-IR spectroscopy also like mixture of these molec- 
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ular forms. However this model can not be used for some quantitative predictions, since 
dimerization described in [22] is considered as a result of dipole-dipole interaction of the 
polar end groups of the monomers. 

Pershin et al. [23] propose a model which takes into account the influence of the molecular 
structure on the macroscopic properties of the nematics, consisting of flexible molecules like 
those HOBA (heptyloxybenzoic acids), OOBA (octyloxybenzoic acids) and NOBA (nony- 
loxybenzoic acids), and predicts isostructural transitions in the nematic phase. These tran- 
sitions already were detected by DSC and impedance spectroscopy [24,25]. On the other 
hand, this theory also can not explain the scalar order parameter vs temperature and dimer, 
monomer concentration temperature variations, because the theory considers the molecu- 
lar conformational order-disorder in a general way and the structural element used as a 
conformational unit is not determined as a concrete part of the molecule. 

Because of these insufficiencies, Barbero et al. presented a theoretical model [1] and it 
is an improved one of the classical mean field theories. It takes into account the internal 
degree of freedom of the molecules due to the variation between the three molecular forms: 
closed dimers, open dimers and monomers. In their model, the energy of the hydrogen 
bonding is one of the quantities connected with the influence of the dimerization process 
on the macroscopic properties. This energy was measured in [16]. It is important to note 
that Barbero et al. have not taken into account the short range forces and fluctuations, 
which leads to a smectic short range order in a coexistence with the nematic long range 
order. Our study can be considered to succeed this disadvantages. So this study is a further 
development of the classical theories, and of the theory proposed by Barbero et al. 

We use the same notation as Barbero et al.; N is the total number of the closed dimers, 
at T = 0; N*, Nm = 2iV*, N^. and are the total number of dissociated dimers, monomers, 
closed dimers and open dimers at a given temperature T. > is the energy necessary to 
break a hydrogen bond of the closed dimer, the activation energy of the chemical reaction 
1 closed dimer— >1 open dimer is then E\ the activation energy of the chemical reactions 1 
closed dimer^2 monomers is 2E and 1 open dimer^2 monomers is E. 
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At a given temperature, the equilibrium distributions of open dimers, closed dimers and 
monomers are given by: 

Uo = exp(-(/i + E)/kBT) 
He = exp{-iJ,/kBT) 
n^ = 2exp(-(// + 2E)//cBT) 

with the condition no + nc + rim/^ — 1, where /i is the chemical potential of the mixture, 

no = No/N, Tic = Nc/N, Um = Nm/N. 

Now we consider the nematic order in the mixture, which is formed by closed dimers, open 
dimers and monomers. Barbero et al. assumed that monomers did not contribute to the 
nematic order. This points out to suppose that the monomers are nearly of spherical shape. 
Although this assumption is not too far from the reality, it could affect the predictions of the 
model proposed by Barbero et al. As can be further, our results is in very good agreement 
with the experimental data. In their study, the interactions among the molecules of the 
closed and open dimers are responsible for the nematic order. They supposed also that the 
closed and open dimers were formed by rod-like molecules. Let n be the mean nematic order 
of the mixture. The angle between the molecular long axis of the closed dimer and n will 
be indicated 9c, and the angle between that of the open dimer and n will be indicated 9o, 
Therefore the total nematic potential is given by 

ViOo, dc) = E J = "^00 + Kc + Kc + (8) 

hj 

in the mean field approximation, where Vpo and Vcc are the mean field potentials acting on 
a molecule of open dimer due to the other open dimers and closed dimers, respectively. Vcc 
and Vco have similar meanings; Vcc and are the potentials acting on a molecule of closed 
dimer due to the other closed dimers and open dimers, respectively. The partial nematic 
mean field potentials Vij are given by [26] 

v,j = -a,,P2{e,) {p,{e,)) (9) 
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in the Maier-Saupe approximation, where the couphng constants ccij depend on the distance 
between the centers of mass of the molecules and on a molecular property. 
Eq.(8) can be rewritten by using Eq.(9) as the following: 

v{eo,e,) = Vo{eo) + v{e,) (lo) 

where 

ViiOi) = - [auSi + a,jS^] P^iOi). (11) 

In Eq.(ll), we use the relation Si = {P2{0i)) for the scalar order parameters of closed {Sc) 
and open {So) dimers, respectively, where ^2(6*1) is the second Legendre polynomial. Then 
the partition function can be given by 

/= dicosOo) [' d{cose,)eM-{Vo + Vc)/kBT) (12) 
Jo Jo 

and the free energy can be written as 

F = U-TSe 

= -NksT ln/+ ItV^q;. .5,5., (13) 

where Sg is the entropy of the system. The scalar order parameters are determined from the 
following relation: 

fi 

where i stands for o, c; open dimers, closed dimers, P = l/(fcsT). Eq.(14) is a self-consistent 
equation and can be solved numerically, and also depends on; a) the temperature T, b) the 
activation energy E, and c) the quantities a^j which contains the strengths of the interaction 
closed dimer-closed dimer, closed dimer-open dimer and open dimer-open dimer. Barbero 
et al. made an analysis and assumed that the concentration dependence of these quantities 
was determined by 



We use the same assumption in this study. The nematic order parameters of the closed (Sc) 
and open (So) dimers can be determined from Eq.(14). If So and Sc are known, then the 
nematic order parameter of the mixture is calculated by using the following relation: 

ric + rio 

Let us now discuss how the order parameter of the mixture is calculated. The partition 
function of the mixture is written as 

"1 rl 



fq = I I d{coseo) d{cosec) exp^(-(K + Vc)/kBT 

= [' /'d(cos^„)rf(cos^e) [l + (l-g)(K + K)ABT]T^ (17) 



within TT. The order parameters are evaluated from Eq.(7), and this equation is written in 
TT formalism as: 

_ /o(dcos^^) pjSi _ 

where i — o,c. If the partial order parameters are evaluated from this equation, the nematic 
scalar order parameter of the mixture is calculated from Eq.(16). The critical value of the 
order parameter, in equilibrium, is determined from Eq.(13). 



II. RESULTS AND DISCUSSION 

The nematic-isotropic transition being first-order, occurs when the orientational 
Helmholtz free energy vanishes, provided the volume change is negleted, and so to investigate 
this transition the free energy is calculated from the interaction potential energy. 

In the Fig.(l), the nematic order parameter of the mixture S vs temperature is shown, 
g = 1 curve represents the Barbero et al. model. It is clear that the curve for q = 1.013 
which is plotted by assuming aij oc , Ucc — 1 eV [26], Uoo — Uco — Uoc — Ucc, and 
E — A.Skcal / mole [16] is in very good agreement with the experimental data. Fig. (2) shows 
the free energy vs the order parameter of the mixture at the phase transition temperature. 
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The curves have two minima; one at s = and the other at s — Sc- s — corresponds to 
the hquid phase, and s 7^ corresponds to the nematic phase. The minimum corresponding 
to the nematic phase for g = 1 occurs at Sc = 0.429. Since The Barbero et al. model is 
based on the Maier-Saupe theory, this is not a surprising result. It is well known that Maier- 
Saupe theory gives a universal value, 0.429, of the order parameter at the nematic-isotropic 
transition temperature. However, the minimum corresponding to the nematic phase for 
q — 1.013 occurs at Sc = 0.442. As can be seen from these results, the minimum of the 
free energy, corresponding to the nematic phase, changes to a higher value as q is increased. 
The experimental results indicate that the critical value of the scalar order parameter in 
nematic-isotropic transition has a value between 0.25 — 0.5 [27]. Although the theory based 
on the Maier-Saupe theory does predict the variation of the order parameter, the quantitive 
predictions of the theory which yields an interaction potential with just one parameter can 
not account precisely for such properties as the temperature dependence of the orientational 
order (Si) [28]. Then generalized form of the standard theory could give a different value of 
the order parameter from 0.429, e.g. 0.442. When this generalized form of the theory based 
on the Maier-Saupe theory is used (e.g. inserting another parameter, q (entropic index) 
coming from the TT), the critical value of the order parameter appears to assume a value 
lying in this range and the temperature dependence of the order parameter is explained 
in very good agreement with the experimental results. As a consequence, with a small 
departure from the standard theory, one could explain the behaviour of the order parameter 
vs temperature and the nematic-isotropic phase transition. 

As becomes clear from the results obtained by employing TT, the present study can also 
be applied to other liquid crystal systems. In this context, a theoretical model, which is 
proposed to investigate the influence of the dimerization process on the nematic ordering, 
has been employed and the theoretical results have been improved by using TT. 
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Figure Caption 



Figure 1. Orientational order parameter as a function of temperature in the mixture for 
q — 1.013 and q — l. Filled squares represent the experimental data[16]. 



Figure 2. The Helmholtz free energy as a function of generalized order parameter in the 
mixture for q — 1.013 and q — l- 
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